Antiproton Production in p + A Collisions at AGS Energies 
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Inclusive and semi-inclusive measurements are presented for antiproton (p) production in proton- 
nucleus collisions at the AGS. The inclusive yields per event increase strongly with increasing beam 
energy and decrease slightly with increasing target mass. The p yield in 17.5 GeV/c p+Au collisions 
decreases with grey track multiplicity, N g , for N g > 0, consistent with annihilation within the target 
nucleus. The relationship between N g and the number of scatterings of the proton in the nucleus is 
used to estimate the p annihilation cross section in the nuclear medium. The resulting cross section 
is at least a factor of five smaller than the free p — p annihilation cross section when assuming a 

.— ». ' small or negligible formation time. Only with a long formation time can the data be described with 

|^_ | the free p — p annihilation cross section. 
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Antiprotons are a topic of great interest in Relativis- 
tic Heavy Ion Physics |l}-|5| because enhanced produc- 
tion may result from the formation of the Quark Gluon 
Plasma (QGP) M and because antiprotons may provide 
an experimental measure of the baryon density of mat- 
ter produced in A + A collisions due to the large p — p 
annihilation cross section [J7| . Comparisons of antiproton 
production in E802 || and E878 || with cascade mod- 
els show that, in fact, more antiprotons are produced 
in these collisions than would be expected from ordi- 
nary hadronic production and the effects of final-state 
absorption. However, the same models have shown that 
these larger p yields could result from increased produc- 
tion and/or decreased absorption. One proposed mech- 
anism for increased production of p is the multiple scat- 
tering of the incident baryons llQ ] . Several mechanisms 
for decreased absorption have been proposed including a 
finite p formation time, "shielding" of the absorption pro- 
cess O] (ARC), and a time delay due to the formation 
of a p — p quasi-bound state |l2|]l3| ] (RQMD). Proton- 
nucleus collisions provide a valuable tool for disentan- 
gling these competing effects and elucidating the dynam- 
ics of antiproton production and absorption in a nuclear 
environment because the density of absorbers (nucleons) 
in the nucleus is well understood. In particular, at en- 
ergies close to the antiproton production threshold, it is 
expected that most antiprotons will be produced from 
the first scattering of the proton in the nucleus. The 
remaining thickness of the nucleus may then simply act 
as an absorber. Previous measurements of p production 
in p + A collisions have been hindered by poor statistics 
and the fact that most p + A measurements are inclu- 
sive. Antiproton yields presented here are shown first 
inclusively, comparing different beam energies and tar- 
gets, and then as a function of the number of scatterings 
of the projectile proton in the nucleus. We focus on our 
semi-inclusive measurement to address questions about 
the first-collision model and in-medium modifications to 
the annihilation cross section. 

E910 is a TPC based p + A experiment with down- 
stream tracking, Cerenkov (CKOV), and time-of- flight 
(TOF) detectors. The E910 apparatus has been de- 
scribed in detail elsewhere |14| , but here we will again 
give a brief description of the time-of-flight (TOF) de- 
tector, used to identify antiprotons, and the trigger used 
for this data set. The TOF is located approximately 
8 m from the target with an active area of approximately 
5x2 m 2 . It is made up of 32 scintillating slats with read- 
out at the top and bottom of each slat. Protons can 
be separated from pions and kaons by more than 2 a 
up to a momentum of 3.0 GeV/c and by more than 1 a 
up to 3.5 GeV/c, where the timing resolution for pro- 
tons is 164 ps. The data presented here were collected 
with a scintillating fiber trigger. The scintillating fiber 
was placed approximately 2 cm downstream of the tar- 
get. A minimum bias interaction is defined as an event 
having two hits in each of the two layers of the scintil- 
lating fiber, and a central interaction is an event with 



a total of 20 or more hits in the fiber. We include only 
central triggers for which there were greater than 2 "grey 
tracks" in the event, where a grey track is defined as a 
"slow" proton or deuteron. A slow proton has momen- 
tum 0.25 < p < 1.2 GeV/c and a slow deuteron has 
momentum 0.5 < p < 2.4 GeV/c. 

For track quality, we require at least 10 hits in the time 
projection chamber (TPC) and that the track originated 
from the event vertex. Although our primary particle 
identification for antiprotons comes from time of flight, 
we additionally use Cerenkov information and ionization 
energy loss in the TPC to reduce background. We require 
that the ionization energy loss is within 3 a of the proton 
dE/dx. In the relativistic rise region, where the pion 
band separates from the proton, the measured dE/dx 
must be greater than 1.5 a from the pion dE/dx. We 
also apply a cut on the Cerenkov ADC which, on average, 
corresponds to requiring less than 0.35 photoelectrons. 
The effect of these cuts on the background can be seen in 
Fig. [l], which will be described later in the text. Quality 
cuts on the hits on the time-of-flight wall (TOF) include 
a cut on the difference in horizontal position between a 
projected track and the center of a hit TOF slat, and 
a cut on the energy deposited on the TOF slat. We 
further require the projected track to have at least 5 hits 
in the 3 drift chambers located between the TPC and the 
TOF wall. We determine the efficiencies of these cuts 
and correct for them in y and px- Tracks are matched 
to the TOF wall with an 90±5% efficiency. A single 
correction factor is applied uniformly over y and px- The 
5% systematic uncertainty in this correction is included 
in the overall errors of the results. 

Figure [fl shows the momentum dependence for nega- 
tive tracks of the difference between the measured time 
of flight and the expected time of flight assuming the 
mass of a proton. The cuts on this distribution are mo- 
mentum dependent and range from 800 ps for low mo- 
mentum tracks to 200 ps for tracks with momenta be- 
tween 3 and 3.5 GeV/c. A momentum-dependent back- 
ground of the identified antiprotons is calculated and sub- 
tracted, amounting to an overall correction of approxi- 
mately 5%. The data have been acceptance corrected 
within our y — px coverage and corrected for the effi- 
ciencies of our cuts as mentioned above. Our coverage 
ranges from 10-800 MeV/c in px and 1-2 in units of ra- 
pidity. The acceptance is largest at low px near rapidi- 
ties between 1.6 and 2.0. The acceptance is limited in 
the low y, high px region by the spatial coverage of the 
TOF wall, while the high y region is limited by the upper 
momentum cut of 3.5 GeV/c. Since a sample of identi- 
fied protons, using loose TOF criteria, is relatively clean 
compared to a sample of loosely identified antiprotons, 
the efficiency of the cuts can be estimated by applying 
the cuts to the sample of identified protons. The over- 
all calculated efficiency ranges from 60-80%, decreasing 
with increasing px- The cut on the difference in hori- 
zontal position of the projected track from the TOF hit 
is most inefficient near the edges of the TOF wall, thus 



affecting primarily the low y, high pt region. The par- 
ticle identification cuts (on dE/dx and Cerenkov ADC), 
which reduce the background, affect the efficiency pri- 
marily for higher momentum antiprotons (high y, high 
Pt region). We neglect the inefficiency due to multiple 
hits in one TOF slat because the slat occupancy is only 
approximately 3%. Even in central events, where the 
multiplicities are larger, the slat occupancy is less than 
3% because the particles tend to shift back in rapidity 
and thus out of the acceptance of the TOF wall. The data 
have also been corrected for trigger biases which are de- 
termined by examining beam-triggered events (unbiased 
events) in two dimensions, the total number of charged 
particles in the event and the number of grey tracks N g in 
the event. Particularly in "central" interaction triggered 
events, but also in "minimum bias" interaction triggered 
events, there is a bias against events with small num- 
bers of charged particles and small N g . This bias can 
be determined from beam-triggered events, by compar- 
ing the distribution of the number of charged particles 
and N g for those beam-triggered events that also passed 
the conditions for an interaction trigger to the unbiased 
distribution. We estimate feeddown from antilambdas, 
by applying the p selection cuts to a set of antilambdas 
identified in the TPC, to be less than 5%. Final (raw) 
event statistics for each data set are shown in Table Q. 
The data sets analyzed include 17.5 GeV/c momentum 
p + Au collisions, 12.3 GeV/c p + Au, 12.3 GeV/c p + Cu, 
and 12.3 GeV/c p + Be. The target thicknesses are 4.5%, 
3.1%, and 2.0% of the interaction length for Be, Cu, and 
Au, respectively. 

Our measure of centrality is defined by the number of 
projectile collisions v, which is derived from the num- 
ber of grey tracks N g in an event. Slow protons and 
deuterons are identified by their measured ionization en- 
ergy loss in the TPC. The momentum cut on the pro- 
tons is 0.25 < p < 1.2 GeV/c and on the deuterons is 
0.5 < p < 2.4 GeV/c. For a class of events with a given 
number of grey tracks N g , we derive the mean number 
of collisions. The details of our method to determine 
< v{N g ) > are described elsewhere fl4J1 . 

Antiproton yields for all 4 data sets are shown in Fig. 
as a function of rapidity and transverse mass. The beam 
energy dependence is seen by comparing the two p + Au 
data sets. The yields increase approximately by a factor 
of 3 from 12.3 to 17.5 GeV/c. The yields for the dif- 
ferent target sizes can also be compared in this figure. 
The transverse mass wit distribution for each target and 
beam momentum is fit to the following exponential, 
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where Co and T are fit parameters. The results of the fits 
are tabulated in Table [fj]. Due to limited statistics, the 
errors on the fit parameters are large. However, compar- 
ing the inverse slope parameter for the largest data set, 
17.5 GeV/c p + Au, with that of the 12.3 GeV/c p + Be 



data set shows a significant increase from the smaller 
target to the larger target. Although the comparison is 
between two data sets with different beam energies, such 
behavior is consistent with more reabsorption in a larger 
target. The beam energy dependence of < pt >, as pre- 
viously measured in p + p collisions |15[] , cannot solely 
account for the observed difference in slope. 

The yields tend to increase with decreasing target size. 
This trend is more evident when the yields are integrated 
over the entire range of measured y (1-2) and pr (10- 
800 MeV/c). These yields, dn/dy with y=l-2 and inte- 
grated over pt, are shown in Fig. ||. There is a 34 ± 22% 
decrease in dn/dy from the Be target to the Au target. 
The 17.5 GeV/c p + Au yield is shown in the same figure 
for comparison. The yield from p + Au at 17.5 GeV/c is 
3.1 times the yield at 12.3 GeV/c. Although the likeli- 
hood of producing antiprotons may be greater in a larger 
nucleus [16| , the likelihood of reabsorption is also greater 
due to the presence of more baryons. At these beam en- 
ergies, we find the effect of increased reabsorption in the 
larger nucleus to be greater than any possible increase in 
production. 

The dependence of the yields on the beam momentum 
can be described by the available kinetic energy squared, 
ft was shown [[t7J that antiproton yields for p + p colli- 
sions at energies near the production threshold can be 
described by 



{KEf = (V^-4m) 2 , 
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where m is the mass of the antiproton. In the reference, 
it is also shown (using phase-space arguments) that this 
dependence can be explained by production through a 
three-body process rather than a four-body process, in- 
dicating the possibility of an intermediate state. Figure [| 
shows that this dependence also describes well our p+Au 
yields. Having established this dependence, we can com- 
pare to data at a different beam energy. 

Figure o shows a comparison between our measure- 
ment of dn/dy and the measurement by E802 fllq] . For 
the purpose of direct comparison, we have restricted our 
y range from 1 to 1.6 and scaled the yields using (KE) 2 
to correspond to the 14.6 GeV/c beam momentum of 
E802. Although the measurements are consistent with 
each other, E802 concluded no target dependence. With 
increased statistics and a larger range in y (Fig. 0), we 
conclude that there is indeed a modest target dependence 
at AGS energies. 

In addition to our inclusive measurements for different 
targets and beam energies, we can use the dependence 
of the p multiplicity on centrality to help disentangle the 
mechanisms of production and reabsorption in the nu- 
cleus. The centrality (or v) dependence of p yields is 
shown in Fig. @. The antiproton yields are measured for 
each value of N g and then plotted versus the mean v for 
a given N g . Since the N g — bin corresponds to a mean 
of approximately two projectile collisions for p + Au, we 
have plotted the first N g bin for p + Be on the same fig- 
ure (scaled from 12.3 to 17.5 GeV/c beam momentum 



using the (KE) 2 relationship discussed above) simply as 
a reference of what the production may be in only one 
p + N collision. Although we may be somewhat biased 
against antiproton production in the N g — bin (the 
antiproton may be preferentially produced together with 
a slow proton due to baryon number conservation), it is 
possible that the increase from N g = to N g — 1 is due 
to contributions to production beyond the first p + N 
collision. However, with the exception of the yield in the 
first p + Au N g bin to the second, we see the mean an- 
tiproton multiplicity decrease as the number of projectile 
collisions increases. This relationship gives insight to the 
amount of nuclear material traversed by the antiproton 
before it is reabsorbed. With a few phenomenological 
assumptions, we present a quantitative measure of the 
survival probability of an antiproton as a function of the 
amount of nuclear material through which it passes. Be- 
cause the beam energy is close to production threshold 
and the antiprotons are strongly peaked at forward an- 
gles, we assume that only the first collisions contribute to 
the production of antiprotons, which are then assumed 
to follow the path of the projectile through the nucleus. 
Since we have conjectured possible contributions to pro- 
duction beyond the first collision, we will discuss the ef- 
fect of such a modification to our assumptions later in 
the text. With this picture of p production, we quantify 
the reabsorption with the following equation, 



a(pA — > pX) — o(pp — > pX)e 
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where <j(pp —* pX) is the antiproton production cross 
section for p + p, a a bs is the "effective" antiproton ab- 
sorption cross section, and a P N is the proton-nucleon in- 
teraction cross section. Since the "v" plotted on the x- 
axis of Fig. ra is simply an average value, v(N g ) : and each 
value of N g actually has a distribution of v values associ- 
ated with it, Pn„{v), we fold the above exponential with 
Pn (v). We determine a a bs by fitting the ^-dependent 
antiproton yields with the following function, 



o(pA -* pX) = a(pp -> pX)P Ng (v)e 
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Folding the distribution, Pjy {y), for discretized values 
of N g results in a step-like behavior of the fit. We show 
smoothed fit functions in the figure. The results of the fits 
with various sets of assumptions are shown in Table III . 
We have done the fits with and without allowing for a 
formation time, 
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during which the antiproton cannot annihilate. In this 
case, the exponent in Eq. becomes — ^hs.( v _ v . — ]_) 

for v > Vf orm + 1, and there is no absorption (a a b s = 0) 
for v < Vform + 1- One should note that the linear rela- 
tionship between Vf orm and Tf orm is not as straightfor- 
ward for very large values of vf orm which rely on fluctu- 
ations in the nuclear density distribution. The quantities 



that are used to calculate a a b s and r/ orm from the fit pa- 
rameters are a mean free path A of 2 fm, a proton-nucleon 
interaction cross section <j p n of 30 mb, and a free annihi- 
lation cross section a ann of 38 mb (at the mean measured 
momentum of 2.5 GeV/c for the antiprotons we detect). 
Using a momentum of 2.5 GeV/c, we calculate 7 and v. 
In the first 3 fits shown in the table, we include the first 
p + Au data point (N g = 0) in the fit, and in the next 
set of 3 fits, we do not include this point (because of the 
initial increase in yield from N g = to N g = 1). In ad- 
dition to removing a possible bias in the first data point 
from the fit, this also allows for production beyond the 
first collision up to the value of v = 2.4 where the fit 
begins. 

The first row shown in Table |lj is the result of a fit 
assuming that the formation time is negligible. The ex- 
tracted <j a bs is significantly reduced relative to a an n (al- 
most by a factor of 10). This fit is shown in Fig. |as a 
dashed curve. In the following 2 fits, we investigate the 
effect of a formation time Tf orm on this result. With no 
constraints on the fit, Tform is very large. However, with 
excessively large Tf orm we lose the ability to uniquely de- 
termine a a bs from the fit. Thus, we constrain Tform and 
Cabs separately in 2 fits. Typical values used in transport 
models for the formation time are 1-2 fm/c [|L9| . Con- 
straining Tf orm to such values results again in a reduced 
absorption cross section. The fit parameters shown in 
fit 2 in the table are for a constraint of Tf orm — 1 fm/c, 
with which one obtains a a a b s — 4.6 ± 0.9 mb. (This fit 
looks similar to the fit with no formation time and is, 
therefore, not shown in the figure.) The other possibility 
is to constrain the absorption cross section to be equal 
to the free annihilation cross section. Such a constraint 
leads to Vf orm = 6.7 ± 0.7, which corresponds to a long 
formation time of 4.9 ± 0.5 fm/c or a formation length 
of approximately 13 fm in the nuclear rest frame. Again, 
such large values of Vf orm rely on density fluctuations 
and thus do not have such a well-defined relationship to 
Tform- The fit is shown by the dotted curve in Fig. |6j. 
The next 3 rows repeat the 3 fits described above, ex- 
cluding the first p + Au data point from the fit. The 
results are qualitatively similar to those when including 
the first data point. Fit 4 in the table, which is shown 
with a solid curve in the figure, again assumes no forma- 
tion time and results in a reduced absorption cross sec- 
tion (approximately 5 times smaller than a ann ). When 
including a formation time of 1 fm/c, fit 5 in the table, 
a a bs is still reduced by approximately a factor of 5. Fi- 
nally, fit 6 shows that constraining a a b s — cr ann results in 
a long formation time, even when excluding the first data 
point. The large discrepancy between a a b s and cr ann , as 
derived from our simple model when the formation time 
is negligible or small (1 fm/c), suggests that the "effec- 
tive" annihilation cross section is very different from the 
free annihilation cross section due to in-medium effects. 
On the other hand, the data can be described by the free 
p — p annihilation cross section and a very long forma- 
tion time. A possible explanation for such a result is the 



formation of an intermediate state where Tf orm can be 
interpreted as the mean lifetime of the state. With such 
a large Tf orm , the antiproton is born only in the late 
stage of the propagation of this state through the nu- 
cleus, leaving little opportunity to get reabsorbed. This 
hypothesis could be tested by measuring the absorption 
of other antibaryons which could proceed through the 
same intermediate state. 

In conclusion, we find the yields dramatically increase 
with increasing beam energy and can be described by 
a dependence on the available kinetic energy squared. 
The observed energy dependence can be understood if 
the antiproton is produced through the decay of an in- 
termediate state |17J . We find a moderate decrease with 
increasing target mass, 34 ± 22% from Be to Au for 
beam momentum 12.3 GeV/c. Finally, we have quan- 
tified the survival probability of an antiproton in the nu- 
clear medium as a function of the number of collisions. 
With this relationship and the assumption of a negli- 
gible or small formation time of 1 fm/c, we find that 
the annihilation cross section is greatly modified within 
the nuclear medium, and that the "effective" absorption 
cross section is a small fraction of the free annihilation 
cross section. On the other hand, a full p — p annihi- 
lation cross section would require a much longer forma- 
tion time than normally assumed. Previous attempts to 
explain a suppression of the annihilation of antiprotons 
within the nucleus include an increased hadron forma- 
tion time ||,|2(J , the formation of a p — p molecule [ |12|Jl3| 
with a finite lifetime, and a "shielding" effect due to the 
presence of mesons pi. All of such phenomenological 
arguments manifest themselves as a delay time during 
which the antibaryon cannot annihilate with a baryon. 
Shielding, in particular, is dependent on the density of 
the nuclear medium and would probably not be a large ef- 
fect in proton-nucleus collisions. Our ^-dependent yields, 
however, show that the reabsorption of antiprotons is al- 
ready greatly suppressed in p + A collisions. Production 
through an intermediate state, which does not get ab- 
sorbed like an antiproton, could also explain a suppres- 
sion of the annihilation of antiprotons. In conclusion, we 
observe a deviation from the expectations of the naive 
first-collision model in which the antiproton is produced 
on-shell in the first collision with a small or negligible for- 
mation time and then interacts with nucleons according 
to the free p — p annihilation cross section. 
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TABLE I. Final (raw) event statistics. 



Target 


P 

(GeV/c) 


No. min. 
bias triggers 


No. central 
triggers 


No. p 
min. bias 


No. p 
central 


Au 


17.5 


2.66 M 


1.06 M 


346 


93 


Au 


12.3 


1.69 M 


0.46 M 


73 


6 


Cu 


12.3 


1.26 M 





84 





Be 


12.3 


1.41 M 





102 






TABLE II. Fit parameters of exponential fits to transverse 
mass distributions. 



Target 


P (GeV/c) 


Co (GeV~V) 


T (MeV/c) 


Au 


17.5 


4.90 ± 0.62 x 10~ 4 


157 ± 34 


Au 


12.3 


1.85 ± 0.91 x 10" 4 


98 ±58 


Cu 


12.3 


2.84 ± 0.95 x 10 -4 


108 ± 47 


Be 


12.3 


4.10 ± 0.92 x 10 -4 


86 ±19 



TABLE III. Fit parameters of antiproton absorption fits. 
Fits 1-3 include the first p + Au data point in the fit, and 
fits 4-6 do not. 



Fit 


(Tabs{mb) 


Tform(fm/c) 


Constraint 


X'/NDF 


1 
2 
3 


4.0 ± 1.6 

4.2 ± 1.6 

38 




1 

4.9 ± 0.5 


Tform U 
Tform -L 
0~abs =38 


6.994/7 
6.562/7 
3.149/7 


4 
5 
6 


6.9 ± 2.2 

6.9 ± 2.4 

38 



1 

4.7 ± 0.7 


Tform — U 
Tform -L 
Cabs =38 


2.005/6 
1.995/6 
1.928/6 
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FIG. 1. Total momenta for negative tracks vs. difference 
between measured flight time and expected flight time in ns. 
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FIG. 3. Target mass dependence of E910 p yields summed 
over 1 < y < 2, T - yields for 12.3 GeV/c beam momentum, 
□ - 17.5 GeV/c. 
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FIG. 2. Target and beam momentum dependence of p + A 
p spectra, a) dn/dy distributions, b) mr distributions, • - 
p+Be, ▼ - p+Cu, ■ - p+Au at 12.3 GeV/c, □ - 17.5 GeV/c 
p+Au. The points from different data sets are offset relative 
to each other in order to distinguish the error bars. 






CO 

0.7 


' 7 — 


1.0 XlO' 4 (KE) 2 


0.6 


: 




0.5 


: 


J 


0.4 


: 


/ 


0.3 


: 


/ 


0.2 


: 


/ 


0.1 




^ttT 


, , i , , , i , , , i , , , i , , , i 


( 


) 8 


10 12 14 16 18 



p (GeV/c) 

FIG. 4. Energy dependence of E910 p + Au p yields, • - p 
yields vs. beam momentum, o - p production threshold. The 
curve shows yields follow KE 2 dependence. 
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FIG. 5. Comparison between E910 12.3 GeV/cp yields (T) 
extrapolated to 14.6 GeV/c and E802 measurements (□). See 
text for details. The points from the 2 experiments are offset 
relative to each other in order to distinguish the error bars. 
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FIG. 6. Dependence of 17.5 GeV/c p+Au p yields on 
(v(N gIcy )) (□). Lines show results of absorption fits. Fit 1 
is shown by the dashed curve, fit 3 the dotted curve, and 
fit 4 the solid curve (see Table |l| for details). Also shown 
is 12.3 GeV/c p+Be data point extrapolated to 17.5 GeV/c 
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